The very-high energy emission from pulsars: a case for inverse 

Compton scattering 



Maxim Lyutikov 
Department of Physics, Purdue University, 
525 Northwestern Avenue, West Lafayette, IN 47907-2036 

Nepomuk Otte 

Santa Cruz Institute for Particle Physics and Department of Physics, 
University of California, Santa Cruz, CA 95060 

Andrew McCann 

Department of Physics, McGill University, Montreal, QC, Canada H3A 2T8 



The observations of gamma-ray emission from pulsars with the Fermi-LAT de- 
tector and the detection of the Crab pulsar with the VERITAS array of Cherenkov 
telescopes at energies above 100 GeV make it unlikely that curvature radiation is 
the main source of photons above GeV energies in the Crab and many other pul- 
sars. We outline a model in which the broad UV-X-ray component and the very 
high energy 7-ray emission of pulsars are explained within the Synchrotron-Self- 
Compton (SSC) framework. We argue that the bulk of the observed radiation is 
generated by the secondary plasma, which is produced in cascades in the outer 
gaps of the magnetosphere. We find that the inverse-Compton (IC) scattering 
occurs in the Klein-Nishina regime, which favors synchrotron photons in the UV 
band as target field for the scattering process. The primary beam is accelerated 
in a modest electric field, with a field strength that is of the order of a few percent 
of the magnetic field near the light cylinder. Overall, in the Klein-Nishina regime 
of the IC scattering the particle distribution in the gap does not evolve towards 
a stationary distribution and thus is intrinsically time-dependent. We point out 
that in a radiation reaction-limited regime of particle acceleration the gamma-ray 
luminosity L 7 scales linearly with the pulsar spin-down power E, L 1 oc E, and 



not proportional to V E as expected from potential-limited acceleration. 
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Introduction 



The recent launch of the Fermi Gamma-Ray Space Telescope and subsequent detection 



of a large number of pulsars (Abdo & et al. 2010) revolutionized our picture of the non- 



thermal emission from pulsars in the gamma-ray band from 100 MeV up to about 10 GeV. 
At even higher energies, in the very-high energy (VHE) band, the detection of the Crab 



pulsar at 25 GeV by the Magic Collaboration (Aliu & MAGIC Collaboration 2008) and 



recently at 120 GeV by the VERITAS Collaboration QAliu fe VERITAS CollaborationpOlTj ) 
in the very-high energy (VHE) band allow to stringently constrain the very-high-energy 
emission mechanisms in the case of the Crab pulsar. In this paper we show that it is very 
difficult to invoke curvature radiation as the dominant radiation mechanism to explain the 
observed emission above 100 GeV and, furthermore, demonstrate that inverse-Compton (IC) 
upscattering of UV photons into the VHE band can explain the observations in the gamma- 
ray band. 



Cheng et al. (1986) were amongst the first to discuss high energy emission from the 
magnetosphere of pulsars. They proposed the outer gap as the location where charged 
particles accelerate to relativistic energies and radiate in the gamma ray band. The outer 
gap model is currently one of the most favored models to explain non-thermal radiation from 
pulsars, [j Based on the idea of the outer gap, geometrical models are very successful in 
explaining the basic features of the observed 7-ray light curves (e.g., Romani & Yadigaroglu 



1995] |Harding et aL]|2008t |Bai fe Spitkovik^|[20ro| . While there seems broad consensus that 
the particle accelerator is located in the outer magnetosphere, the radiation physics remain 
controversial. One of the preferred radiation mechanisms, which is believed to dominate 



the observed gamma-ray emission, is curvature radiation Romani (1996) (see also Cheng 



et al. (1986, 2000); Takata et al. (2008); Tang et al. (2008)). Possible importance of the 



IC scattering was mentioned previously (e.g., Romani 1996), but was never considered the 
primary emission mechanism for the very high energy photons. 

This paper is structured in the following way. In $2] we demonstrate that the recent 



results obtained with Fermi (Abdo feet al. 2010), Magic (Aliu fe MAGIC Collaboration 2008 ) 



and especially VERITAS QAliu fe VERITAS Collaboration|201l| ) make it highly unlikely that 
curvature emission is the main radiation mechanism of photons above 10 GeV energies from 
the Crab pulsar. In §|3] we show that inverse-Compton scattering by secondary particles in 
the outer gaps is broadly consistent with the observed luminosity in the very-high-energy 
band. 



1 Below, for the order-of-magnitude estimates, by "outer gaps models" we imply generic "outer magneto- 
sphere models" . 
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Limits on curvature radiation 



2.1. Crab pulsar 



Curvature radiation is a widely discussed process to explain the observed gamma-ray 
emission from the magnetosphere of pulsars. In this section we discuss the difficulty of invok- 
ing curvature radiation as the emission process that explains the observed pulsed emission 
from the Crab pulsar above 100 GeV. 

We assume an outer gap scenario with the accelerating electric field being parallel to the 



magnetic field Cheng et al. ( 1986 ). In the electric field a beam of charged particles accelerates 



- hereafter primary beam - that has a particle density, which is of the order of the Goldreich- 
Julian density uqj ( Goldreich fc Julian||1969 ). The primary beam loses a significant amount 
of its energy through various radiative processes of which the curvature emission and the 



IC- induced pair production are the dominant ones Arons (1983); Cheng & Ruderman (1977). 



The pair-production process results in the formation of a second population of particles - 
hereafter secondary plasma -, which has a higher particle density than the primary beam but 
a smaller bulk Lorentz factor. Within this outer-gap framework we derive a general upper 
limit of the break in the curvature radiation spectrum that is emitted by particles within the 
outer gap of the Crab pulsar. The limit we obtain is independent of the particular details 
of the acceleration mechanism of the primary beam. In our argument we follow a similar 
approach that has been applied before in the discussion of the synchrotron emission from 



pulsar wind nebulae by de Jager et al. (1996) and Lyutikov (2010). 



Within the outer gap, the charged particles follow the curved magnetic field lines and, 
therefore, emit curvature-radiation photons. The curvature radiation spectrum emitted by 



monoenergetic particles has a break at energy ( Zheleznyakov 1996) 

3 



f 6r 
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where R c is the curvature radius of the magnetic field lines, and 75 is the Lorentz factor of 
the radiating particles. 

An upper limit of 7& is set by the constraint that while the particles accelerate they 
radiate and, therefore, the maximum value of 7j is obtained when acceleration gains are 
balanced by radiative losses, i.e. the radiation reaction limit. Under the assumption that 
the accelerating electric field E is a fraction 77 < 1 of the magnetic field B, the acceleration 
gain is ecrjB, with e is the electron charge. The radiation reaction limit is then reached if: 
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where the losses due to curvature radiation are given on the right side. Using Eq. ^ it 
follows from Eq. Q that 

'3\ 7/4 /— ( 5\ 3/4 

\-^\ hcyRc ( v— ) (3) 



*) 

The radius of curvature R c can be expressed in units of the light cylinder Rl, R c = ^Rl = 
£cP/(27r), where P is the period of the pulsar and £ is a dimensionless scaling parameter. 
If, furthermore, B is replaced by the radial distribution of the magnetic field of a dipole 
B = B NS (R NS /R) 3 Q , where B NS is the magnetic field on the surface of the neutrons star 
and Rns the stars surface, then it follows that: 
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= 150 GeV r/ 3/4 = 5 GeV ^ 
9 x 10 7 r ? 1 / 4 v ^ = 3 x 10y.2 4 Vf- ( 4 ) 



On the very right side we parametrized r/ = 10 2 ?7-2, assuming that the electric field is 



a few percent of the magnetic field. In the current outer-gap models (Cheng et al. 1986 



2000 



Takata et al. 



2008 



Tang et al. 2008) an electric field of E\\ « (QBr 2 )/(cR c ) ~ 0.15 



is predicted, while in the models of Cheng et al. (2000); Takata et al. ( |2008 ); Tang et al. 



(2008) the accelerating field is one order of magnitude smaller. 



In the radiation reaction limited regime, the maximum energies of photons emitted by 
curvature radiation is determined by the maximum energies of the electrons. The result is 
a break in the spectrum and an exponentially falling flux above the break. The gamma-ray 



spectrum of the Crab pulsar has a break at about 6 GeV (Abdo & et al. 2010), which is 



formally consistent with the break we predict for an electric field that is a few percent, i.e. 
?7_2 ~ 1- However, the observed gamma-ray flux above the break is not exponentially falling, 
which is expected if the break would be due to curvature radiation produced by the electrons 
with the highest energies. The non-exponential cutoff can be explained if the electric field 
is larger, i.e. rj ~ 1, or if a different emission mechanism dominates above the break energy. 



2.2. Other pulsars 

In the gamma-ray band a few dozen pulsars have been detected With the Fermi-LAT 
detector. The spectral energy distributions of these pulsars in the 7-ray band are very similar 



2 Very close to the light cylinder the toroidal magnetic field induced by the poloidal currents becomes 
important; we neglect here these effects. 
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and can be characterized by a flat spectral component between 100 MeV and a few GeV 
(Thompson et al. 1999 Abdo &: et al.poTO ) and a spectral break in the GeV region. In Fig. 
[I] we compare the observed spectral breaks of these pulsars with the predicted ones from Eq. 
(J4J) by calculating the ratio of the observed Eb r and predicted spectral break e^r for each of 
the 46 pulsars reported in the first Fermi catalogue ( Abdo fc et al.||2010 ). For the calculation 
of the break energy e^r we used the extreme case of i] = ( = 1. If the spectral break is due 
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Fig. 1. — Ratio of the observed break energies Eb r for 46 pulsars to the maximum predicted 
for curvature radiation eb r , which is given by Eq. Q with r\ = £ = 1 

to curvature radiation and the electric field in the gap is much less than one, 77 << 1, as it 
is expected in present outer gap models, the ratio should be much smaller than one. This is 
indeed the case for the majority of the pulsars, including the Crab pulsar. However, for a 
significant number of pulsars the ratio is close to one and for one pulsar, PSR J1836 + 5925, 
the ratio is even larger than one. In order to explain the spectral break for these pulsars as 
a result of curvature radiation an accelerating electric fields is required that is close to or 
even larger than magnetic fields . 



Two possible interpretations of these results are: (i) the observed spectral break is due 
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to curvature radiation by the electrons with the highest energies. For the Crab pulsar a 
new component dominates above the break and explains the non-exponential cutoff. In this 
interpretation, the pulsars, for which the ratio E^je^. is close to unity, can be explained 
by statistical outliers (uncertainties on are not taken into account in Fig. [I]), or that a 
different emission mechanism dominates at high energies that influences the measurement of 
Ef, r . (ii) The gamma-ray emission above ~ GeV energies is due to one single emission process, 
which is not curvature radiation. In this case the spectral break reflects the underlying 
particle distribution. 



3. IC model of the high-energy gamma-ray emission from the Crab pulsar 

3.1. Outline of the model 

In this section we outline the key features of an SSC model that is able to explain the 
high energy emission of the Crab pulsar. Observationally, the spectral energy distribution 
(SED) of the Crab pulsar has a broad peak in the 10-100 keV range with a luminosity 



Lx ~ 10 36 ergs 1 (cf. Fig. 9 in Kuiper et al. 2001) which is a few percent of the pulsar's 
spin-down power of ~ 5 x 10 38 ergs _1 . Between 10 MeV and the spectral break at a few 



GeV the SED is flat and has a luminosity L 7 of a few xl0 34 ergs 1 (Abdo 2010). Above the 
spectral break at > 150 GeV the luminosity is ~ 10 33 erg s ~ 1 . 

We identify the broad soft UV- X-ray peak in the SED of the Crab pulsar as a syn- 
chrotron (or possibly cyclotron) emission from the secondary plasma boosted by the large 
parallel velocities of emitting particles. This creates target photons for IC scattering both by 
the primary beam and by the secondary plasma. As we demonstrate below, the IC scattering 
by the secondary plasma is broadly consistent with the observations. 



3.2. IC scattering by the primary beam 

In this section we discuss the inverse Compton scattering by the primary beam in the 
outer gap. We use, like in the previous section, E = i]B = 10~ 2 ?7_2-£> for the accelerating 
electric field and t]gR z lc with t]g = O-lrjc-i for the effective emitting volume. In order to 
simplify our calculations, we separate the broad UV-X-ray peak into two component: A low 
energy component that covers the UV band with a luminosity of Luv ~ 10 34 ergs _1 Z34 and 
typical photon energies of e so ft = leVe^o, and a high energy component that is centered 
around the X-ray peak with a luminosity of Lx ~ 10 36 ergs _1 /36 and typical photon energies 
of e S oft — lkeVex,3- The need to separate the broad-band component into two comes from 
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the strong dependence of the IC scattering in the KN regime on the energy of the photon . 

The properties of Inverse Compton scattering strongly depend on whether the scattering 
occurs in the Thompson regime or in the Klein-Nishina (KN) regime. In which regime the 
scattering takes place is determined by the Lorentz factor of the scattering particle and the 



energy of the upscattered photon (Blumenthal & Gould 1970). For a given photon energy e 



the scattering takes place in the KN regime if the Lorentz factor ^kn is larger than 

Ikn = ]— « 1.2 x 10 5 e^ « 1.2 x 10 2 e^ 3 (5) 

These are fairly modest Lorentz factors considering the above estimate of the maximum 
Lorentz factor that can be achieved in the outer gap in the radiation reaction limit, Eq. Q. 
It can, therefore, be concluded, that inverse Compton scattering by the primary beam takes 
place in the KN regime. 

Adding losses due to inverse Compton scattering in the extreme KN limit into the 
balanced gain loss equation ^ results in a net energy loss of (Blumenthal & Gould 1970 



Schlickeiser & Ruppel 2010) 



4 = ec " B - Tr<* tk) " s O u -*™ (6) 

where U so f t is the energy density of the target photon field, e so f t is the typical energy of a soft 
photon, and <7j is the Thompson cross-section. Note that both the acceleration term and the 
decelerating IC term are independent of the energy of the particle. Thus, if curvature losses 
were negligible, particles are either accelerated or decelerated without reaching a steady 
solution. Only in the presence of curvature radiation is it possible to achieve a steady-state 
particle distribution. 

In order to better understand how curvature radiation and inverse-Compton scattering 
contribute to the radiation loss of the primary beam in Eq. (JBJ we compare the two. In this 
comparison we assume the curvature radiation-limited Lorentz factor Eq. Q, and justify our 
choice post factum by showing that curvature radiation and IC losses in the Crab pulsar are 
about equal. The soft photon luminosity that results in IC losses in the KN regime which 
are similar to curvature radiation losses is 

_ BnsKs^uv _ f 10 35 ergs- 1 4 V0 7/_ 2 

-ksoft.crit — V ~2 - 1 ln 41 -1 2 K' ) 

e A \ 10 ergs e xz r\^2 

The minimum luminosity of the target photon field in the UV that is needed to achieve IC 
losses similar to curvature radiation losses Eq. ([7]) is about the same as the observed UV 
luminosity. The upscattering of soft X-ray photons, even though the X-ray flux is higher 
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than the UV flux, does not contribute much to the radiative loss of the primary beam (KN 
suppression) because the observed X-ray luminosity is five orders of magnitude below the 
critical luminosity. 

The conclusion that the IC upscattering of UV photons and curvature radiation con- 
tribute about equal to the total loss of the primary beam means that both processes also 
contribute equally to the emitted power in the gamma-ray band. However, the two processes 
produce very different spectral features. As we have shown before, curvature radiation pho- 
tons can only be emitted with energies up to a few GeV for reasonable electric fields and 
curvature radii, see Eq. Q. The spectrum of the IC upscattered photons, on the other 
hand, extends to much higher energies. This can be shown by assuming again that curva- 
ture radiation and IC losses are about equal, in which case the maximum Lorentz factor can 
still be estimated with Eq. Q. The maximum energy of the upscattered photons, e 7 , is 
then given by the maximum electron energy: 

7/4 pl/4 p 3 / 4 

e 7 - 7b m e c 2 = (37r)V4^V /4 V^^F = 15TeV^ 4 V! (8) 

While the maximum photon energy produced by IC scattering depends on the maximum 
electron energy, the total power emitted by IC scattering is independent of the electron 
energy. Instead the total power is determined by the low-energy target photons field. Due 
to the steeply falling IC cross-section in the KN regime with increasing energy of the target 
photons oc e~ 2 , the maximum IC power L^n might not be determined by the peak luminosity 
in the spectral energy distribution of the target photons but be at lower energies: 



J KN,b 



2\ 2 

U so{t a T c x n GJ x r] G R\ c (9) 

Csoft / 



Application to the Crab pulsar yields that the primary beam produces an IC luminosity by 
upscattering the X-ray photons with keV energies, luminosity of Lx ~ 10 36 erg s _1 , that is, 

L KN , x = 5xlQ 29 r lG! _ 1 e x 2 3 (10) 

This is much lower then the IC luminosity produced by upscattering the UV photons with 
eV energies, Ljjv ~ 10 34 erg s^ 1 : 

L K n,uv = 5 x 10 33 T] G _i . (11) 

The above is an estimate of the peak power. The average luminosity is lower by at least 
one order of magnitude. Thus, we conclude that the IC scattering by the primary beam is 
unlikely to be the origin of the VERITAS signal. 



-9- 



3.3. Gamma-ray emission from the secondary plasma 



In the previous section we discussed the gamma-ray emission produced by the primary 
beam. In this section we discuss the gamma-ray emission by the particles that are produced 
in pair cascades of the particles in the primary beam, the secondary plasma. 

We recall that the primary beam has a density uqj and a Lorentz factor 7& Q. As 
nomenclature for the secondary plasma we use n p for its density and 7 P for its Lorentz factor. 
We assume energy equipartition between the primary beam and the secondary plasma (the 
assumption of equipartition between the primary beam and secondary plasma is justified 
in the polar cap models; Daugherty fc Harding 1996, we assume a similar parametrization 



here). From equipartition it follows that n P 7 P = ncj'jb- The two particle populations are 
connected through the pair cascading process, i.e. n p = X p Ugj, where A = IOOA2 is the 



multiplicity factor of the secondary particles. Multiplicities of the order A 
in outer gap models (e.g. 



10 are typical 



A ~ 10 4 - 10 b (Takata et al. 2010). 



Wang & Hirotani 2011), but can also reach much higher values, 



In our picture of a radiation-reaction limited acceleration of the primary beam it follows 
that the Lorentz factor of the secondary plasma is given by 



7 P «7 6 /A = 3xl0 5 ? 7 i / 2 yiA 2 - 1 



(12) 



This Lorentz factor is above the minimum Jkn- Therefore, IC scattering by the secondary 
plasma takes place in the KN regime and we can use the same relations that we have derived 
in the previous section for the emission produced by the primary beam. (For multiplicities 
mush higher than the assumed A ~ 100 the scattering by UV photons occurs in the Thompson 
regime. Overall, the convolution of the electron and the soft photon spectrum requires 
detailed radiative calculations which include global magnetospheric models and anisotropic 
angular distributions of the photons.) The maximum energy of IC photons produced by the 
secondary plasma is (cf. Eq. Q) 



,c 2 = 150GeVr ] 1/ ^X 2 1 



and the peak luminosity of the IC scattered UV photons is (cf. Eq. (11) is 



Lrn, p — XL 



KN,b 



4 x 10 35 r? G ,-i%vo A 2 



(13) 



(14) 



Both the energy (13) and the peak luminosity (14) are consistent with the VERITAS de- 



tection. Thus, IC up-scattering of UV photons by the secondary plasma can explain the 
observed pulsed emission from the Crab pulsar above 100 GeV. We leave a more detailed 
calculation of the spectrum to a future paper (we expect that the overall spectrum will 
depend on the distributions both in parallel m and perpendicular p± momenta). 
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The secondary plasma also produces synchrotron photons with energies that can be 
estimated, e.g., using Doppler-boosted cyclotron emission: 

e X , P = hu Blp = ^ /4 V^ A ^f 7/4 = SkeVr^v^ 1 . (15) 

This roughly coincides with the energy where the Crab pulsar emits most of its power (in 
fact, the Crab emits most of it's power around 100 keV). 

To produce the observed synchrotron luminosity L s N p (e 2 /c)u^ , ~fj_'-fp ~ 10 36 ergs _1 
(N p is the total number of secondary particles in the magnetosphere, 7j_ is a typical transverse 
Lorentz factor) one requires 

N « 5 x 10 32 e^^ o7 - 2 2 

7± = ^= (16) 

This demonstrates that for the chosen parameters 7_i_ ~ 1, the soft emission occurs in the 
cyclotron regime, and it also shows that the overdensity 

A = f f m C ) = 60 ^ e uv,o, (17) 

n GJ\^LC) 

is consistent with our assumption of A 2 ~ 1. 

We, therefore, conclude that emission from the secondary plasma is not only able to 
explain the observed gamma-ray emission above 100 GeV by upscattering UV photons but 
it also explains the bulk of the X-ray emission. An obvious modification is required to 
this simplified picture to include the relativistic momenta of the secondary particles that is 
transverse to the magnetic field lines and results in synchrotron and not cyclotron emission 
as we assumed. 



4. Expected X-ray-7-ray correlations 

Within the framework of the SSC model the power emitted by IC is related to the power 
of the seed photons. Photons of different energies that are emitted by the same particles 
should in principle produce similar pulse profiles. In our model one expects, therefore, that 
the pulse profiles in X-ray and in gamma-rays are similar because the secondary plasma 
emits synchrotron radiation in X-rays and IC scatters UV photons into the VHE band. And 
indeed, the ratio of the amplitudes of the two pulses in the pulse profile of the Crab pulsar 
changes consistently in the X-rays / soft gamma-ray band and in the high energy gamma-ray 
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band. In X-rays the main pulse dominates over the inter pulse. The ratio changes towards 
higher energies and reverses in the soft gamma-ray band at about 1 MeV. Similarly, the main 



pulse dominates at 100 MeV (see, e.g., also for the pulse profiles at lower energies Abdo & 



et al. 2010) while at 120 GeV the inter pulse clearly dominates over the main pulse Aliu & 



VERITAS Collaboration (2011). 



In addition, as we have argued that though IC losses may be energetically dominant (or 
similar to curvature emission), in the KN regime they do not lead to a equilibrium distribu- 
tion of Lorentz factors. Hence we expect highly non-stationary magnetospheric plasma flows. 
This will lead to highly non-stationary radiative properties. Since within the SSC model the 
soft and hard photon fields are related, we might expect some 7-ray - X-ray correlation. 
Though it is the soft UV photons that are scattered to the GeV energies, and, formally, one 
expects UV-GeV correlation, since X-rays and UV form a continuous spectral distribution, 
one also expects X-ray - GeV correlation as well. Thus we expect short time-scale statistical 
correlation between X-ray and 7 rays photons. 



5. Dependence of the 7-ray luminosity on the spin-down power in the 

radiation-reaction-limited regime 

Here we discuss the dependence of the 7-ray luminosity on the spin-down power in the 
radiation-reaction limit. Generalizing Eq. (|2]), the total luminosity radiated by a primary 
beam of Goldreich- Julian density n G j = BQ/(2irec) in the radiation reaction limited regime 
is 

L c = ecr]Bn GJ r] G R 3 LC (18) 

where rjcB? LC is the volume occupied by the radiating particles. Replacing B with a dipole 
field Bns * (Rns/R) 3 at the light cylinder R = Rlc, it follows that 

Lc w mG B% s & N j& w (ig) 

where Esd ~ - MS ^^—, is the pulsar spin-down power. 

Thus, in the radiation-reaction-limited regime, the gamma-ray luminosity is proportional 
to the spin-down power, L c oc E S d- This differs from the commonly used L c oc a/ E S d scaling, 
which results if the maximum particle energy is not limited by radiation reaction but by the 
electric potential and most of the energy is radiated away once the particle is outside of the 
accelerating region. This is the case in polar cap models. In these models a beam with a 
particle density equal to the Goldreich- Julian density loses energy N oc n GJ r 2 PC oc 



SD 



see e.g., Zhang & Harding 2000), where rp G is the radius of the polar cap. The same 
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square-root scaling has been extended to outer gaps, assuming that the emitting volume is 



proportional to the volume within the light cylinder radius (Hirotani et al. 2003). 



The expected linear proportionality (19) of the 7-ray luminosity is valid in the radiation 
reaction limit, i.e. if the dominant radiation processes depend on the particle energy. This 
is the case, for example, for curvature radiation or inverse-Compton (IC) scattering in the 
Thompson regime. And it is not the case for IC scattering in the Klein-Nishina regime, 
where the radiative losses are independent of the particle energy (see Eq. [5]), and, therefore, 
the acceleration is not limited by radiation losses. However, as we argued in £j3j there are 
good reasons to believe that particle acceleration is indeed limited by radiation reaction. 

We note that testing our prediction is complicated by the large uncertainty of the geo- 
metrical parameter rjc, the effective emission volume, which depends on the pulsar inclination 
angle, the angle between the rotation axis, and the line of sight. It may also depend on the 
period of the pulsar through the microphysics of the acceleration. 



6. Discussion 



The recent detection of the Crab pulsar above 100 GeV by VERITAS (Aliu & VERITAS 



Collaboration 2011) changes our picture of high energy gamma-ray emission from pulsars. 
Even though the breaks in the energy spectra of most pulsars are consistent with curvature 
radiation in the radiation reaction limit. For some pulsars exceptional conditions on the 
accelerating electric fields are required to explain the observed cutoff energies with curvature 
radiation. In particular, the pulsed emission from the Crab pulsar above 100 GeV can 
only originate from curvature radiation if extreme assumptions are being made about the 
pulsar's magnetosphere. The observation that the flux above 100 GeV is in agreement with 
an extrapolation of the flux from the GeV regime argues in favor of one emission mechanism 
being dominant below and above the spectral break. 

Thus, there are two somewhat independent arguments against curvature radiation as 
the dominant source of GeV photons: (i) in many pulsars the observed break energy is too 
high; (ii) the Crab pulsar energy spectrum above the break is inconsistent with what is 
expected if the break is due to curvature radiation from particles in the radiation reaction 
limited regime. 

More precise measurements of the energy spectra of all 7-ray pulsars, but especially for 
those with low break energies, may be decisive for further progress: for those pulsars for 
which the break energy is consistent with curvature radiation and moderate electric fields; 
it is then expected that the energy spectrum above the break follows an exponential cut-off. 
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If, however, the energy spectra above the breaks are better described by power laws like for 
the Crab, it argues agains curvature radiation. 

In this paper we demonstrated that inverse Compton scattering in the Klein-Nishina 
regime by the secondary particles results in an overall consistent picture with observations. 
The key features of our model are (i) A population of primaries that is accelerated in a 
modest electric field, which is a fraction rj of the magnetic field strength near the light 
cylinder with a typical value of rj is 10~ 2 . The suppression of the scattering cross-section in 
the Klein-Nishina regime (and the corresponding lower radiation loss rate of electrons) allows 
primary leptons to be accelerated to very high energies with hard spectra, (ii) The gain in 
energy of the primaries in the electric field is balanced by similar curvature radiation and 
IC losses (radiation reaction limited); (iii) The secondary plasma is less energetic, but more 
dense and has approximately the same energy content as the primary beam. The secondary 
plasma is responsible for the soft UV-X-ray emission via synchrotron/cyclotron emission and 
the high energy 7-ray emission that extends to hundreds of GeV via the inverse Compton 
process. The IC emission from the primary beam extends well into the TeV regime but will 
be difficult to detected due to the low predicted fluxes. 

Finally, we argued that in the radiation reaction-limited regime the 7-ray luminosity of 



pulsars should scale linearly with the spin-down energy, Eq. (19). The coefficients of this 
proportionality depend both on the overall geometry of the magnetosphere (e.g., inclination 
angle of the magnetic dipole with respect to the axis) through the parameters rjc and on the 
electric field in the gap through the parameter 77 (which, in turn, depends on microphysics 
of the acceleration precesses) . 

This prediction is in contrast to the currently assumed scaling of the 7-ray luminosity 
with the available potential, oc a/ Esd- Observationally, when compared with the scaling 
of oc a/ Esd, all models underpredict the luminosity of pulsars and thus fail to describe the 
observed population (e.g., Pierbattista et al.j|2011 ) (for alternative interpretation of data see 



Watters fc Romam]|2011 ). The proposed linear scaling of the 7-ray luminosity with the spin 



down energy naturally predicts more energetic pulsars. 

Here we outlined a framework to explain the non-thermal radiation from gaps in the 
magnetosphere of pulsars. More detailed calculations of the emitted energy spectra are 
needed. A major complication in including the IC loses in the radiation codes results from 
the fact that in the KN regime, an accelerating electric field of a given strength does not 
lead to a fixed energy of a particle. This means that a particle is either accelerated or 
decelerated depending on the photon density and the value of the electric field and does not 
reach a steady energy. This implies that in this regime acceleration is highly non-stationary. 
The addition of curvature radiation can, however, establish a steady state. Thus, curvature 
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radiation, even if not dominating the total gamma-ray luminosity, may dictate the particle's 
final energy. 

A number of additional factors must be taken into account to construct a comprehen- 
sive model of the higher energy emission. Most important is the intrinsically non-isotropic 
distribution of soft photons. A more detailed structure of the magnetic field lines within the 
magnetosphere need to be taken into account, including modifications due to magnetospheric 
currents. Also, particle trajectories may not exactly follow the magnetic field lines due to 
various drift effects. An important modification could be the IC scattering of the surface 



thermal emission closer to the surface of the neutron star in the slot gaps (Arons 1983) (in 



comparison, Crab does not show any thermal component Weisskopf et al. 2004) 



Our model is based on the assumption that emission is generated within the light cylin- 
der. The main argument for this is that Fermi pulsar profiles are well fitted with geometric 
models and that the profile of the very high energy emission is correlated with the lower 
energies. This disfavors the models that advocate the emission from the wind zone (e.g., 



Bogovalov & Aharonian 2000 Kirk et al. 2002). 
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